maThis is REGCOIL_PM,
a regularized least-squares method for designing permanent magnets for stellarators.
 Successfully read parameters from regcoil_nml namelist in regcoil_in.benchmark.
 Resolution parameters:
   ntheta_plasma  =   96
   ntheta_coil    =   96
   nzeta_plasma   =    4
   nzeta_coil     =    4
   mpol_magnetization =   24
   ntor_magnetization =    0
   ns_magnetization =    1
   ns_integration   =    2
 Not assuming stellarator symmetry
 We will use the following values of the regularization weight lambda:
 0.000E+00 1.000E-34 4.833E-34 2.336E-33 1.129E-32 5.456E-32 2.637E-31 1.274E-30 6.158E-30 2.976E-29 1.438E-28 6.952E-28 3.360E-27 1.624E-26 7.848E-26 3.793E-25 1.833E-24 8.859E-24 4.281E-23 2.069E-22 1.000E-21 4.833E-21 2.336E-20 1.129E-19 5.456E-19 2.637E-18 1.274E-17 6.158E-17 2.976E-16 1.438E-15 6.952E-15 3.360E-14 1.624E-13 7.848E-13 3.793E-12 1.833E-11 8.859E-11 4.281E-10 2.069E-09 1.000E-08
 Initializing plasma surface.
 Reading FOCUS format data from file:plasma
 Number of modes for Bnormal read from FOCUS file:           1
 Plasma surface area: 8.290E+02 m^2. Volume: 2.900E+02 m^3.
 G =    0.0000000000000000       ; curpol =    1.0000000000000000     
 Done initializing plasma surface. Took    2.30000000E-02  sec.
 Initializing coil surface.
   Building a plain circular torus.
   Evaluating coil surface & derivatives:   8.00000038E-03  sec.
 s_integration:  0.21132486540518713       0.78867513459481287     
 s_weights:  0.50000000000000011       0.50000000000000011     
 s_magnetization  0.50000000000000000     
 s_magnetization_weights   1.0000000000000000     
 interpolate_magnetization_to_integration:
   1.0000000000000000     
   1.0000000000000000     
 Maximum allowed uniform d before singularity:  2.9000E+01
 Coil surface area: 1.184E+03 m^2. Volume: 5.918E+02 m^3.
 Done initializing coil surface. Took    9.99999978E-03  sec.
 Using Bn coefficients in FOCUS file.
 Initializing basis functions
 Done. Took   0.00000000     sec.
 Building inductance matrix.
   Number of OpenMP threads:           4
 Done. Took   1.59099996     sec.
 Computing g and RHS_B:
  js=   1 j_RZetaZ= 1 dgemm took  1.000E-03 sec, rest took  0.000E+00 sec.
  js=   1 j_RZetaZ= 2 dgemm took  1.000E-03 sec, rest took  0.000E+00 sec.
  js=   1 j_RZetaZ= 3 dgemm took  0.000E+00 sec, rest took  0.000E+00 sec.
 Total time to compute g and RHS_B:   2.00000009E-03 sec.
 About to form matrix_B.
  js=   1, ks=   1, j_RZetaZ= 1, k_RZetaZ= 1 dgemm took  0.000E+00 sec.
  js=   1, ks=   1, j_RZetaZ= 2, k_RZetaZ= 1 dgemm took  0.000E+00 sec.
  js=   1, ks=   1, j_RZetaZ= 2, k_RZetaZ= 2 dgemm took  0.000E+00 sec.
  js=   1, ks=   1, j_RZetaZ= 3, k_RZetaZ= 1 dgemm took  0.000E+00 sec.
  js=   1, ks=   1, j_RZetaZ= 3, k_RZetaZ= 2 dgemm took  0.000E+00 sec.
  js=   1, ks=   1, j_RZetaZ= 3, k_RZetaZ= 3 dgemm took  0.000E+00 sec.
 Total time for matrix_B:   0.00000000     sec.
 Regularization term:   1.00000005E-03 sec.
 Optimal LWORK:        9408
 Solving system for lambda= 0.000E+00 (  1 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  2.58720278802226E+27, method 2 =  2.58720278802226E+27, relative difference = -2.12490422642228E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.75056430748701E-15, method 2 =  8.60779891809216E-10, relative difference = -9.99993609154804E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 2.751E-15,  chi2_M: 2.587E+27
   max(B_n): 3.341E-09,  max(M): 4.233E+15,  min(M): 1.261E+13
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 1.000E-34 (  2 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  4.97103553216427E+19, method 2 =  4.97103553216427E+19, relative difference = -4.11986594492989E-16
   2 methods of computing chi2_B that should agree: method 1 =  6.31783845742711E-17, method 2 =  1.61435309564695E-11, relative difference = -9.99992172946861E-01
   Diagnostics:    0.00000000      sec.
   chi2_B: 6.318E-17,  chi2_M: 4.971E+19
   max(B_n): 4.515E-10,  max(M): 4.599E+11,  min(M): 1.035E+10
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 4.833E-34 (  3 of  40)
   Additions:    1.00000005E-03  sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  4.87678283987662E+19, method 2 =  4.87678283987662E+19, relative difference =  0.00000000000000E+00
   2 methods of computing chi2_B that should agree: method 1 =  6.28946818488622E-17, method 2 =  2.48405740421731E-11, relative difference = -9.99994936145799E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 6.289E-17,  chi2_M: 4.877E+19
   max(B_n): 4.491E-10,  max(M): 4.557E+11,  min(M): 1.075E+10
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 2.336E-33 (  4 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  4.87143104989705E+19, method 2 =  4.87143104989704E+19, relative difference =  5.88574480605776E-16
   2 methods of computing chi2_B that should agree: method 1 =  6.30629707630821E-17, method 2 =  6.63931132294238E-11, relative difference = -9.99998100317705E-01
   Diagnostics:    0.00000000      sec.
   chi2_B: 6.306E-17,  chi2_M: 4.871E+19
   max(B_n): 4.498E-10,  max(M): 4.555E+11,  min(M): 1.077E+10
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 1.129E-32 (  5 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  4.85088786860040E+19, method 2 =  4.85088786860040E+19, relative difference = -2.53314451557203E-16
   2 methods of computing chi2_B that should agree: method 1 =  4.95459820866363E-17, method 2 =  5.27506927028298E-11, relative difference = -9.99998121505785E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 4.955E-17,  chi2_M: 4.851E+19
   max(B_n): 4.062E-10,  max(M): 4.546E+11,  min(M): 1.094E+10
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 5.456E-32 (  6 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  1.44599743213358E+19, method 2 =  1.44599743213358E+19, relative difference = -6.37345530164721E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.41569391159225E-16, method 2 =  1.36424205265939E-12, relative difference = -9.99792478612509E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.416E-16,  chi2_M: 1.446E+19
   max(B_n): 6.082E-10,  max(M): 2.130E+11,  min(M): 3.802E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 2.637E-31 (  7 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  2.41361211779415E+19, method 2 =  2.41361211779415E+19, relative difference =  3.39408305899908E-16
   2 methods of computing chi2_B that should agree: method 1 =  3.11356885888628E-16, method 2 =  2.49542608798947E-11, relative difference = -9.99975046105144E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 3.114E-16,  chi2_M: 2.414E+19
   max(B_n): 8.852E-10,  max(M): 3.511E+11,  min(M): 1.752E+10
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 1.274E-30 (  8 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  1.06621718480683E+16, method 2 =  1.06621718480683E+16, relative difference =  9.37895218956890E-17
   2 methods of computing chi2_B that should agree: method 1 =  8.38499180955897E-19, method 2 =  5.11590769747272E-13, relative difference = -9.99996721997908E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 8.385E-19,  chi2_M: 1.066E+16
   max(B_n): 4.641E-11,  max(M): 4.888E+09,  min(M): 5.197E+08
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 6.158E-30 (  9 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  1.54558345824497E+15, method 2 =  1.54558345824498E+15, relative difference = -9.70507281245921E-16
   2 methods of computing chi2_B that should agree: method 1 =  5.23177391714868E-19, method 2 =  4.54747350886464E-13, relative difference = -9.99997699044145E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 5.232E-19,  chi2_M: 1.546E+15
   max(B_n): 3.670E-11,  max(M): 2.347E+09,  min(M): 5.070E+07
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 2.976E-29 ( 10 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  1.69792445828043E+15, method 2 =  1.69792445828043E+15, relative difference =  0.00000000000000E+00
   2 methods of computing chi2_B that should agree: method 1 =  1.12655104221842E-19, method 2 =  4.54747350886464E-13, relative difference = -9.99999504537735E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.127E-19,  chi2_M: 1.698E+15
   max(B_n): 1.705E-11,  max(M): 2.380E+09,  min(M): 5.414E+07
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 1.438E-28 ( 11 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  1.19304539259105E+16, method 2 =  1.19304539259105E+16, relative difference =  8.38191074882912E-17
   2 methods of computing chi2_B that should agree: method 1 =  2.98343612786165E-20, method 2 = -5.68434188608080E-13, relative difference =  1.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 2.983E-20,  chi2_M: 1.193E+16
   max(B_n): 8.770E-12,  max(M): 4.747E+09,  min(M): 1.946E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 6.952E-28 ( 12 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  3.12077644121229E+14, method 2 =  3.12077644121230E+14, relative difference = -3.00406010382410E-16
   2 methods of computing chi2_B that should agree: method 1 =  3.10912521640934E-22, method 2 =  6.82121026329696E-13, relative difference = -9.99999999088395E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 3.109E-22,  chi2_M: 3.121E+14
   max(B_n): 8.981E-13,  max(M): 7.237E+08,  min(M): 2.821E+08
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 3.360E-27 ( 13 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  2.94456350171360E+14, method 2 =  2.94456350171360E+14, relative difference = -5.30638921215555E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.42290300740238E-23, method 2 =  4.54747350886464E-13, relative difference = -9.99999999937420E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.423E-23,  chi2_M: 2.945E+14
   max(B_n): 2.047E-13,  max(M): 5.757E+08,  min(M): 4.273E+08
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 1.624E-26 ( 14 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  2.93074271287233E+14, method 2 =  2.93074271287233E+14, relative difference =  0.00000000000000E+00
   2 methods of computing chi2_B that should agree: method 1 =  6.93255985095112E-25, method 2 =  5.68434188608080E-13, relative difference = -9.99999999997561E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 6.933E-25,  chi2_M: 2.931E+14
   max(B_n): 6.350E-14,  max(M): 5.384E+08,  min(M): 4.529E+08
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 7.848E-26 ( 15 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  2.93004390380973E+14, method 2 =  2.93004390380973E+14, relative difference =  2.13307383956724E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.34319553042827E-24, method 2 =  6.82121026329696E-13, relative difference = -9.99999999996062E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.343E-24,  chi2_M: 2.930E+14
   max(B_n): 7.538E-14,  max(M): 5.300E+08,  min(M): 4.575E+08
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 3.793E-25 ( 16 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  2.93002508719397E+14, method 2 =  2.93002508719397E+14, relative difference = -5.33271884540885E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.99237268117013E-23, method 2 =  7.95807864051312E-13, relative difference = -9.99999999924797E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 2.992E-23,  chi2_M: 2.930E+14
   max(B_n): 2.836E-13,  max(M): 5.286E+08,  min(M): 4.583E+08
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 1.833E-24 ( 17 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  2.93002401482776E+14, method 2 =  2.93002401482776E+14, relative difference =  3.19963247828572E-16
   2 methods of computing chi2_B that should agree: method 1 =  6.96471751067715E-22, method 2 =  7.95807864051312E-13, relative difference = -9.99999998249648E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 6.965E-22,  chi2_M: 2.930E+14
   max(B_n): 1.328E-12,  max(M): 5.283E+08,  min(M): 4.585E+08
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 8.859E-24 ( 18 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  2.93002393018666E+14, method 2 =  2.93002393018666E+14, relative difference = -2.13308838047675E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.62574692323923E-20, method 2 =  5.68434188608080E-13, relative difference = -9.99999942799117E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.626E-20,  chi2_M: 2.930E+14
   max(B_n): 6.381E-12,  max(M): 5.282E+08,  min(M): 4.586E+08
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 4.281E-23 ( 19 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  2.93002392854898E+14, method 2 =  2.93002392854898E+14, relative difference = -5.33272095417251E-16
   2 methods of computing chi2_B that should agree: method 1 =  3.79693185509630E-19, method 2 =  6.82121026329696E-13, relative difference = -9.99998886728426E-01
   Diagnostics:    0.00000000      sec.
   chi2_B: 3.797E-19,  chi2_M: 2.930E+14
   max(B_n): 3.080E-11,  max(M): 5.282E+08,  min(M): 4.586E+08
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 2.069E-22 ( 20 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  2.93002392778741E+14, method 2 =  2.93002392778741E+14, relative difference = -1.06654419111172E-16
   2 methods of computing chi2_B that should agree: method 1 =  8.86837101404566E-18, method 2 =  6.25277607468888E-13, relative difference = -9.99971634214535E-01
   Diagnostics:    0.00000000      sec.
   chi2_B: 8.868E-18,  chi2_M: 2.930E+14
   max(B_n): 1.488E-10,  max(M): 5.282E+08,  min(M): 4.586E+08
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 1.000E-21 ( 21 of  40)
   Additions:    1.00000005E-03  sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  2.93002392449838E+14, method 2 =  2.93002392449838E+14, relative difference =  4.26617676923577E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.07139181390602E-16, method 2 =  6.82121026329696E-13, relative difference = -9.99392845871354E-01
   Diagnostics:    0.00000000      sec.
   chi2_B: 2.071E-16,  chi2_M: 2.930E+14
   max(B_n): 7.192E-10,  max(M): 5.282E+08,  min(M): 4.586E+08
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 4.833E-21 ( 22 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  2.93002390861919E+14, method 2 =  2.93002390861919E+14, relative difference = -1.06654419808905E-16
   2 methods of computing chi2_B that should agree: method 1 =  4.83818955859720E-15, method 2 =  6.82121026329696E-13, relative difference = -9.85914186907469E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 4.838E-15,  chi2_M: 2.930E+14
   max(B_n): 3.476E-09,  max(M): 5.282E+08,  min(M): 4.586E+08
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 2.336E-20 ( 23 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  2.93002383187717E+14, method 2 =  2.93002383187716E+14, relative difference =  3.19963267807066E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.13006614879194E-13, method 2 =  8.52651282912120E-13, relative difference = -7.65948965699619E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.130E-13,  chi2_M: 2.930E+14
   max(B_n): 1.680E-08,  max(M): 5.282E+08,  min(M): 4.586E+08
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 1.129E-19 ( 24 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  2.93002346098842E+14, method 2 =  2.93002346098842E+14, relative difference = -1.06654436102904E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.63951932709944E-12, method 2 =  3.24007487506606E-12, relative difference = -1.02142346447214E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 2.640E-12,  chi2_M: 2.930E+14
   max(B_n): 8.118E-08,  max(M): 5.282E+08,  min(M): 4.586E+08
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 5.456E-19 ( 25 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  2.93002166851000E+14, method 2 =  2.93002166851000E+14, relative difference = -1.06654501350127E-16
   2 methods of computing chi2_B that should agree: method 1 =  6.16517810944215E-11, method 2 =  6.23003870714456E-11, relative difference = -5.23271183248792E-03
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 6.165E-11,  chi2_M: 2.930E+14
   max(B_n): 3.923E-07,  max(M): 5.282E+08,  min(M): 4.586E+08
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 2.637E-18 ( 26 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  2.93001300560999E+14, method 2 =  2.93001300560999E+14, relative difference =  0.00000000000000E+00
   2 methods of computing chi2_B that should agree: method 1 =  1.44000962712926E-09, method 2 =  1.44063960760832E-09, relative difference = -2.18693922003728E-04
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.440E-09,  chi2_M: 2.930E+14
   max(B_n): 1.896E-06,  max(M): 5.282E+08,  min(M): 4.586E+08
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 1.274E-17 ( 27 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  2.92997113896020E+14, method 2 =  2.92997113896019E+14, relative difference =  7.46594384808962E-16
   2 methods of computing chi2_B that should agree: method 1 =  3.36341333163708E-08, method 2 =  3.36345920004533E-08, relative difference = -6.81868253441377E-06
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 3.363E-08,  chi2_M: 2.930E+14
   max(B_n): 9.164E-06,  max(M): 5.282E+08,  min(M): 4.586E+08
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 6.158E-17 ( 28 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  2.92976881301222E+14, method 2 =  2.92976881301222E+14, relative difference =  8.53309649859247E-16
   2 methods of computing chi2_B that should agree: method 1 =  7.85545414415533E-07, method 2 =  7.85546205861465E-07, relative difference = -5.03755427879950E-07
   Diagnostics:    0.00000000      sec.
   chi2_B: 7.855E-07,  chi2_M: 2.930E+14
   max(B_n): 4.429E-05,  max(M): 5.281E+08,  min(M): 4.586E+08
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 2.976E-16 ( 29 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  2.92879128120189E+14, method 2 =  2.92879128120190E+14, relative difference = -3.20097920946854E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.83420298593408E-05, method 2 =  1.83420306711923E-05, relative difference = -2.21309082947357E-08
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.834E-05,  chi2_M: 2.929E+14
   max(B_n): 2.140E-04,  max(M): 5.280E+08,  min(M): 4.585E+08
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 1.438E-15 ( 30 of  40)
   Additions:    1.00000005E-03  sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  2.92407382818554E+14, method 2 =  2.92407382818553E+14, relative difference =  6.41228679634086E-16
   2 methods of computing chi2_B that should agree: method 1 =  4.27728546286428E-04, method 2 =  4.27728546924300E-04, relative difference = -7.45650455198101E-10
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 4.277E-04,  chi2_M: 2.924E+14
   max(B_n): 1.033E-03,  max(M): 5.276E+08,  min(M): 4.581E+08
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 6.952E-15 ( 31 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  2.90143459935180E+14, method 2 =  2.90143459935180E+14, relative difference =  4.30821359984904E-16
   2 methods of computing chi2_B that should agree: method 1 =  9.91318353665907E-03, method 2 =  9.91318353720771E-03, relative difference = -2.76723782237357E-11
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 9.913E-03,  chi2_M: 2.901E+14
   max(B_n): 4.975E-03,  max(M): 5.255E+08,  min(M): 4.564E+08
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 3.360E-14 ( 32 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  2.79564148669399E+14, method 2 =  2.79564148669399E+14, relative difference =  9.50139713064986E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.23100273762528E-01, method 2 =  2.23100273763407E-01, relative difference = -1.97087982675430E-12
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 2.231E-01,  chi2_M: 2.796E+14
   max(B_n): 2.359E-02,  max(M): 5.157E+08,  min(M): 4.481E+08
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 1.624E-13 ( 33 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  2.35767960814205E+14, method 2 =  2.35767960814205E+14, relative difference =  0.00000000000000E+00
   2 methods of computing chi2_B that should agree: method 1 =  4.39453104057435E+00, method 2 =  4.39453104057486E+00, relative difference = -5.82076636874067E-14
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 4.395E+00,  chi2_M: 2.358E+14
   max(B_n): 1.046E-01,  max(M): 4.729E+08,  min(M): 4.120E+08
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 7.848E-13 ( 34 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  1.21214953581237E+14, method 2 =  1.21214953581237E+14, relative difference =  3.86709713736639E-16
   2 methods of computing chi2_B that should agree: method 1 =  5.27678196001318E+01, method 2 =  5.27678196001322E+01, relative difference = -4.17429125847527E-15
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 5.277E+01,  chi2_M: 1.212E+14
   max(B_n): 3.606E-01,  max(M): 3.379E+08,  min(M): 2.966E+08
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 3.793E-12 ( 35 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  2.16266271140804E+13, method 2 =  2.16266271140804E+13, relative difference =  8.12800114750924E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.19874827734594E+02, method 2 =  2.19874827734595E+02, relative difference = -1.09873664324160E-15
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 2.199E+02,  chi2_M: 2.163E+13
   max(B_n): 7.315E-01,  max(M): 1.436E+08,  min(M): 1.260E+08
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 1.833E-11 ( 36 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  1.50445011552229E+12, method 2 =  1.50445011552229E+12, relative difference =  9.73673859230261E-16
   2 methods of computing chi2_B that should agree: method 1 =  3.57241877391936E+02, method 2 =  3.57241877391936E+02, relative difference = -3.97793641074477E-16
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 3.572E+02,  chi2_M: 1.504E+12
   max(B_n): 9.294E-01,  max(M): 3.800E+07,  min(M): 3.331E+07
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 8.859E-11 ( 37 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  7.24078822077349E+10, method 2 =  7.24078822077350E+10, relative difference = -7.37568343257593E-16
   2 methods of computing chi2_B that should agree: method 1 =  4.01592171527695E+02, method 2 =  4.01592171527695E+02, relative difference = -4.24635410430713E-16
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 4.016E+02,  chi2_M: 7.241E+10
   max(B_n): 9.845E-01,  max(M): 8.344E+06,  min(M): 7.313E+06
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 4.281E-10 ( 38 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  3.17880642691117E+09, method 2 =  3.17880642691117E+09, relative difference =  7.50025472086493E-17
   2 methods of computing chi2_B that should agree: method 1 =  4.11796981058842E+02, method 2 =  4.11796981058842E+02, relative difference =  1.38037483214783E-16
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 4.118E+02,  chi2_M: 3.179E+09
   max(B_n): 9.968E-01,  max(M): 1.749E+06,  min(M): 1.532E+06
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 2.069E-09 ( 39 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  1.36809239241861E+08, method 2 =  1.36809239241861E+08, relative difference =  3.26757782071378E-16
   2 methods of computing chi2_B that should agree: method 1 =  4.13957036830388E+02, method 2 =  4.13957036830388E+02, relative difference = -2.05975791459122E-16
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 4.140E+02,  chi2_M: 1.368E+08
   max(B_n): 9.993E-01,  max(M): 3.628E+05,  min(M): 3.179E+05
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 1.000E-08 ( 40 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  5.86361313895958E+06, method 2 =  5.86361313895957E+06, relative difference =  3.97077089050903E-16
   2 methods of computing chi2_B that should agree: method 1 =  4.14406104285241E+02, method 2 =  4.14406104285242E+02, relative difference = -4.80089371164008E-16
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 4.144E+02,  chi2_M: 5.864E+06
   max(B_n): 9.999E-01,  max(M): 7.511E+04,  min(M): 6.582E+04
 Volume of magnetization region:  1.18E+00 meters^3.
 REGCOIL_PM complete. Total time=   1.6749999523162842      sec.
 You can run regcoilPlot regcoil_out.benchmark.nc to plot results.
